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Abstract 

We study the dynamics of roots of /'(C) t), where /(£, t) is a locally 
univalent polynomial solution of the Polubarinova-Galin equation for 
the evolution of the conformal map onto a Hele-Shaw blob subject to 
injection at one point. We give examples of the sometimes complicated 
motion of roots, but show also that the asymptotic behavior is simple. 
More generally we allow f'((, t) to be a rational function and give 
sharp estimates for the motion of poles and for the decay of the Taylor 
coefficients. We also prove that any global in time locally univalent 
solution actually has to be univalent. 
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1 Introduction 



Polynomial, rational and logarithmic solutions of the Polubarinova-Galin 
equation for the conformal map onto a growing or shrinking Hele-Shaw blob 
of viscous fluid have been studied in many papers the last few decades, see for 



example [19], [2], [8], [13], [20], [1], [7], [14], [9], [18], [22], [11], [12], and also 



[6]. In [19] D. Bensimon and B. Shraiman set up the dynamical equations, 
in the polynomial case, for how the zeros of the derivative of the conformal 
map move in the complex plane, and they proposed studying the dynamics 
of the zeros. As far as we know nobody has yet undertaken this task to any 
substantial extent. The purpose of the present paper is to start investigations 
in such a direction. 

The general setting is that f(-,t) : D — > Q(t) C C is the time dependent 
normalized (f(0,t) = 0, f'(0,t) > 0) conformal map onto the fluid region 
which has a source of strength q(t) > at the origin. The evolution is 
then described by the Polubarinova-Galin equation 



We focus mainly on the case that the derivative g((, t) = /'(£, t) is a rational 
function, a property which is preserved in time under (1.1). Then g((,t), 
and hence f{(,t), is completely determined by the zeros ui(t), . . . , u m (t) and 
poles Ci(^)> • • • Xn(t) of g((,t), both located outside the unit disk, together 
with an overall scale factor b(t): 



The motion of the poles Q(t) has previously been studied, for example 
it is known [20], [6] that they always move away from the origin. In the 
present paper we obtain more precise estimates of their speed and locations. 
The zeros Uk(t) show up a more complicated behavior than the poles, for 
example they may occasionally move towards the origin, and even reach the 
unit circle, with possible break down of the solution as the result. Zeros may 
also collide with each other or escape to infinity in finite time. However the 
asymptotic behavior is simple: if the solution does not break down, then as 
t — > oo the zeros move to infinity arranged asymptotically as the corners of 
a regular polygon of growing size. 



Re f((,t)(f>((,t) =q(t) forCGdD. 



(1.1) 



g((,t) = b(t) 



nr=i(c -"*(*)) 
n;=i(c-ao)' 
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Besides the zeros and poles we also study the asymptotics of the coeffi- 
cients in the power series 

oo 

/(C,f) = I>(0C'' 

3=1 

in the case that g((,t) is rational. The leading coefficient ai(t) > increases 
to infinity as t — > oo, in fact can be normalized for example so that a\(t) = e* 
(which then fixes q(t)), while the others tend to zero very quickly: 

as t — )> oo. 

Much of the behavior of solutions of the Polubarinova-Galin equation 
(1.1) can be understood in terms of harmonic moments and quadrature iden- 
tities. These provide enough conserved quantities to make the Hele-Shaw 
problem fully integrable, hence in principle algebraically solvable. In prac- 
tice it is not that easy because of the degree of complexity of the integrals 
of motion when expressed for example in terms of the Taylor coefficients of 
f((,t). However the conservation laws allow for natural concepts of weak 
solutions, and with such relaxed forms of solutions it is possible to let ze- 
ros of g{(,t) penetrate the unit circle and enter the unit disk. Non locally 
univalent solutions of this kind will be studied in a forthcoming paper [3], 
and part of the aim of the present paper is to set the stage for these future 
investigations. 

The organization of the paper is as follows. Section 2 contains all nec- 
essary preliminaries, in particular we set up the notations to be used. In 
Section 3 we derive the dynamical equations for the zeros and poles of g((, t) 
by first writing (1.1) as an equation in g alone (see (3.3)) and then identify- 
ing the residues in this formula (Theorem 3.1). When everything is spelled 
out (Theorem 3.2 and 3.3) one gets a rather involved system of ordinary 
differential equations for Uk(t), (j(t), b(t). From Theorem 3.3 we obtain 
precise estimates of the speeds of the poles (Corollary 3.1), and a kind of 
conservation law for the dynamics (Proposition 3.2). 

The Polubariova-Galin equation (1.1) is not sensitive for loss of univalence 
of the solution as long as it stays locally univalent. However, loss of univalence 
is always followed by loss also of local univalence and break down of the 
solution at a later, but finite, time (Theorem 4.1). In [3] we will show how it 
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in such cases still can be continued as a kind of weak solution spreading on 
a branched Riemann surface. 

In Section 5 we give several examples (Example 5.1, 5.2, 5.3) of the be- 
havior of roots u>k(t) in the polynomial case, for example collision of roots, 
plus two theorems asserting the previously mentioned asymptotic behavior. 
In the first (Theorem 5.1) it is assumed a priori that the solution is global in 
time, while in the second (Theorem 5.2) only sufficiently strong assumptions 
on the initial data are made. Section 6, finally, contains precise estimates for 
the poles Q(t) in the rational case (Theorem 6.1) as well as the asymptotic 
estimates (1.3) of the Taylor coefficients of f((,t) (Theorem 6.2; see also 
Lemma 5.1). There are also some sharper estimates (e.g., Corollary 6.2) for 
the case that the sequence of harmonic moments contains gaps. 

2 Preparatory material 
2.1 List of notations 

We first list some general notations which will be used in the paper. 
. D = {( G C : |C| < 1}, D(a,r) = {C G C : |C - a\ < r}. 

• dm = dm(z) = dxdy = ^dzdz (z = x + iy), area measure in the z 
plane. 

• u* = =, for u eC. 

. h*(0 = h(l/() = E™ =1 where fc(£ t) = £™ i b^. 

(There is a slight ambiguity in this notation: we have = l/£ if £ is 
considered as a point, whereas /*(C) = 1/C f° r the function f(() = (.) 

• /(C,0 = f/(C,0,/'(C,0 = f/(C,*)- 

• With E C C any set which contains the origin, 

{/ : / is analytic in some neighborhood of E}, 
{/GO(£):/(0) = 1 /'(0)>0} ) 
{/eOoonn(£):/V0on E}, 
{/ G 0\ 0CVL (E) : / is univalent (one-to-one) on E}. 



0(E) = 

rm(-S) 

0\ 0CU (E) = 
O UIliv (E) = 
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If a(t), b(t) are positive functions a(t) ~ b(t) will mean that there exist 
constants < c < C < oo such that 

- bit) - 

for all t. 



2.2 Basic equations 

In the paper we shall study certain aspects of the Polubarinova-Galin 
equation, which describes the evolution of a simply connected Hele-Shaw 
fluid blob in the complex plane driven by injection or suction at one point, 
chosen to be the origin. A smooth map t >->■ f(-,t) G O un i v (JB>) will be called 
a (univalent) solution of the Polubarinova-Galin equation if it satisfies 



Re 



/(C,*)C/'(C0 =?(*) forCG^D (2.1) 



in the pointwise sense. Here q(t) is a real- valued continuous function, which 
is given in advance and which represents the strength of the source/sink at 
the origin. Typically q = ±1, which corresponds to injection (plus sign) or 
suction (minus sign) at a rate 2n. Since the transformation t \-> — t changes 
q to — q in (2.1) it is enough to discuss one of the cases q > and q < 0. In 
general we shall take q > 0. 

Equation (2.1) expresses that the image domains Q(t) = f(D,t) evolve in 
such a way that 

— / hdm = 2irq(t)h(0) (2.2) 



for every function h harmonic in a neighborhood of Q(t). This means that 
the speed of the boundary dQ(t) in the normal direction equals q(t) times 
the normal derivative of the Green's function of Q(t) with a pole at z — 0. 
The terminology Laplacian growth is also used for this kind of evolution. The 
equivalence between (2.1) and (2.2) follows from the general formula 

d 
~dt 



[ V dm= [ <p(f((, t))Re [/(£ t)C/'(C, 01 & (C = e ie ), (2.3) 
7n(t) L J 



valid for any smooth evolution t >->■ /(•, t) G O univ (D) and for any smooth test 
function tp in the complex plane. The derivation of (2.3) is straightforward 
(and omitted). 
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On choosing h(z) = z k , k = 0, 1, 2, . . . in (2.2) it follows that the harmonic 
moments 



M k (t) = - f z k dm = f f((,t) k f*((,t)f((,t)d( (2.4) 
n Jn(t) 27T1 J m 

are conserved quantities, except for the first one, which is related to q by 
|M (t) = 2q(t). Thus 

M (t) = M (0) + 2Q(t), 
where Q(t) is the accumulated source up to time t > 0: 

Q(t) = [ q(s)ds. (2.5) 
Jo 

The preservation of Mi, M 2 , . . . characterizes Laplacian growth in the simply 
connected case. 

One may consider the equation (2.1) on different levels of generality. It is 
natural to keep the normalization /(0) = 0, f'(0) > 0, in fact the coupling to 
(2.2) depends on this, but (2.1) makes sense for any / e O norm (D), at least as 
long one makes sure that q(t) = whenever a zero of /' appears on <9IB>. In a 
forthcoming paper [3] we shall deal with this general case, while in the present 
paper we shall only consider locally univalent functions, / G 0\ ocn (JB>). The 
mathematical treatment of (2.1) in this case is exactly the same as in the 
'physical' case / G (9 univ (D). We shall then speak of a locally univalent 
solution of the Polubarinova-Galin equation. 

When / G O iocu (B), then //£/' G (9(0) and equation (2.1) can be solved 
for / by dividing both sides with |C/'| 2 - The result is an equation which we 
shall refer to as the Lowner-Kufarev equation, namely 

f(C,t) = cf(c,t)P(c,t) (CgD), (2.6) 

where P((,t) is the analytic function in D whose real part has boundary 
value q(t)\f'((, t)\~ 2 and which is normalized by ImP(0, t) = 0. Explicitly 

^(C>*) is g iven b y 
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The right member in (2.7) is the Poisson-Schwarz integral 

V,(C)= r^^MO) (CGD) (2.8) 
Jo e ~~ s 

for the measure 

d ^=Xfwm d ^ (2 - 9) 

If g(t) > then \x is positive and Re P > in D. 

Remark 2.1. In the general literature on univalent functions, "Lowner- 
Kufarev equation" usually refers to any equation of the type (2.6) without any 
particular coupling between / and P. The Hele-Shaw case is characterized 
by the feed-back relation (2.7) between / and P. 

Expanding / in a power series, 

f((,t)=a l (t)( + a 2 (t)( 2 + ..., (2.10) 

it follows from (2.6) that 

di(t) 



ai(t) 



= P(0,t) = J dfl(;t). (2.11) 



In particular, if q(t) is chosen so that /i becomes a probability measure, 
namely 

1 1 r d i (212) 

q(t) 2nJ |/'(e* f)| 2 ' 1 ' ' 

then P(0,t) = 1 and ai(t) = Oi(0)e*. 

When / G Ciocu(D) then P e 0(D), in fact the right member of (2.6) 
extends analytically as far as / does (see [2]). We shall keep the notation 
P = P((, t) also for the analytic extension of the Poisson integral beyond D. 

As a general notation throughout the paper, we set 

g(C,t) = f(C,t). (2.13) 

The function g in fact turns out to be more fundamental than / itself. Of 
course, / can be recaptured from g by 



,t)= [ g((,t)d(. 
Jo 



Most of the paper will deal with the case that g is a rational function, 
or perhaps better to say, g(()d( is a rational differential, in other words an 
Abelian differential on the Riemann sphere. If g has residues then / will 
have logarithmic poles, besides ordinary poles. The terminology Abelian 
domain for the image domain f2 = f(B) has been used [21] for this case. 
Alternatively one may speak of f2 being a quadrature domain (see [5] for 
the terminology), which in the present case means that a finite quadrature 
identity of the kind 

- r . I rij- 1 

/ h(z)dxdy = J2 C J / h ( z ) dz + Yl Yl ^^{zj) (2.14) 

J Q j=l j=0 k=l 

holds for integrable analytic functions h in Q. Here the points Zj G Q are 
fixed (i.e., independent of h), with specifically z = 0, the Cj, are fixed 
coefficients, and the 7j are arcs in Q with end points among the Zj. This 
sort of structure is stable under Hele-Shaw flow because, as is seen from 
(2.2), what happens under the evolution is only that the right member is 
augmented by the term 2irQ(t)h(0). 

When g — f is rational we shall write it on the form 

Here m > n = X^=i n j> \Cj\ > 1 an( i repetitions are allowed among the ujk, 
Q to account for multiple zeros and poles. Then, with the argument of b(t) 
chosen so that g(0,t) > 0, / G Oi ocu (D) if and only if \u k \ > 1, \Q\ > 1 for 
all k and j. The assumption m > n means that g(()d(, as a differential, has 
at least a double pole at infinity, which the Hele-Shaw evolution in any case 
will force it to have because the source/sink at the origin creates a pole of / 
at infinity. 

The form (2.15) is stable in time, with the sole exception that when m = n 
the pole of / may disappear at one moment of time. Then the value of m 
drops below n at this moment. This can easily be explained in terms of the 
quadrature identity (2.14): under the Hele-Shaw evolution there will be one 
time dependent term, namely (a i + 2nQ(t))h(0), and this may vanish for 
one value of t. If there are no terms a fc/i < - A: ~ 1 - ) (0) with k > 2 this will cause 
/ to temporarily lose its pole at infinity. The phenomenon is illustrated in 
Example 3.2 below. 
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The rightmost member of (2.15) will be used when we need to be explicit 
about the orders of the poles. The convention then is that Ci>--->C< are 
distinct and rij > 1. Thus n = Y^j=i n h an d m the f un sequence 0, . . . , C„, 
the tail Cm-i> • • • ? Cn wm be repetitions of (some of) the 0, • • • , Ct according 
to their orders. In equations (2.14) and (2.15), £ and the rij are the same. 

For later use we here also exhibit log g and its derivatives: 

m n 

hgg(Ct) = log 6(*) + $>g(C - w fc (f)) - £>g(C - (2.16) 

fe=i j=i 

i^>-S-S^& + S* (2 - 17) 

9 m 1 ™ 1 

We remark that log \g\ can be viewed as the logarithmic potential of the 
charge distribution which puts positive unit charges at the zeros of g, negative 
unit charges at the poles (added up according to multiplicities). This is the 
charge distribution corresponding to the divisor (g) of g, which is defined as 
the formal linear combination 

m n 

(g) = 1 • M (0) -(m-n). (oo). (2.19) 

k=i j=i 

(This is the divisor of g as a function, the divisor of gd( as a differential has 
two additional poles at infinity.) One main subject of the present paper is 
the study of the dynamics of (g) under Hele-Shaw flow. 

Note that log \g\, and even its restriction to <9B, contains all information 
of /. In fact, given u = \og\g\ on <9D we can extend it harmonically to D, 
then form its harmonic conjugate v, normalized by v(0) = 0, and finally 
define / by 

/(C) = / exp(u(z) + \v(z))dz. 
Jo 

Note also that (log |<?|)|aB is a "free" function, i.e., is subject to no constraints 
besides regularity (real analyticity is needed to start a Hele-Shaw evolution 
in a full neighborhood of t — 0, see e.g. [16]). 
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The structure (2.15) of g means that / is of the form 

e £ rij-l , m-n+1 

/(c, t) = Yl e i lo §(c - cm) +EE + E ^)c fc . (2.20) 

i=i j=± k=i ^ ywj fc=0 

The coefficients ej, which are the residues of g(()d(, will not depend on 
time (this is a consequence of (2.1)). If they are not all zero, then / is 
single- valued only outside a suitable system of 'cuts' in C \ D connecting 
the logarithmic poles. Note that if X^=i e j then / automatically has a 
logarithmic pole at infinity, hence one of the cuts has to reach infinity in this 
case. The relationship between (2.20) and (2.14) for = Sl(t) = f(B,t) is 
that Zj = f((*(t),t), 1 < j < t, z = /(00V) = f(0,t) = 0, and that the 
above mentioned 'cuts' can be taken to be the reflections in the unit circle 
of the arcs /~ 1 (7 J )- Moreover, Y^j=i n j = n an d n o = m — n + 2, hence 
X]j=o n i = 171 + 2. One may also note that b = (m — n + l)d m _ n+1 . The 
two systems of coefficients, {cj, a^} and {ej, Cjk, d^}, are related by nonlinear 
equations (see e.g. [21]). 



3 Dynamics of zeros and poles in the rational 

case 

In this section we set up the dynamical equations for zeros and poles in the 
locally univalent rational case. We take / and g to be of the form (2.20), 
(2.15), so that the divisor of g is given by (2.19). When g is rational also P 
is rational. When we refer to P outside D we shall always mean this rational 
function, or more generally the analytic continuation of P\® (when it exists). 
First we write the dynamical equations in a general form. 

Theorem 3.1. Under the above assumptions, the evolution of the divisor 
(g) under Hele- Shaw flow governed by the L owner- Kufarev equation (2.6) is 
given by 



u k {t) = -Res ?=Wfe 
CM = Res f=(: 



(P((,t) 



g((,t) 
g'(C,t) 
g((,t) 



(3.1) 
(3.2) 
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valid whenever Uk, Cj ^ s a simple zero/pole (respectively) of g. 

If bJk is a multiple zero, with say Ue = Uk for I in an index set K C 
{l,2,,...,m} (containing k), then the individual to kit) need not exist, but 
we have instead 



d 



— = -Res 



g(C,t) 



For multiple poles the Q do exist, and we have 



nj(j(t) = Res, 



C=Ci 



g((,t) 



for l<j<£. 

Proof. From (2.6), i.e., / = C,gP, we get g = (CgP)' and hence 
§- t ^9((,t) = (P((,t)^\ogg((,t) + ^((P((,t)). 



(3.3) 



At first this equality holds in D, but since both members are rational func- 
tions it becomes an identity between two rational functions. 

From (2.17) we see that the Uk and Q are the residues of log g. Since the 
last term in (3.3), being a total derivative, has no residues we immediately 
get the equations for Uk and Q. The statement concerning multiple roots 
and poles also follows easily. □ 

Remark 3.1. It should be mentioned that multiple zeros never survive for 
any period of time, only collisions can occur. Multiple poles are however 
stable, they never split or collide. 

In addition to depending on the zeros and poles, g in (2.15) also depends 
on the factor b. The m + n + 1 complex parameters u±, . . . , u m , Ci, • • • , (n, b 
are connected to the coefficient a±(t) = g(0,t) > by 



l-iy 



n 



fc=i 



CLi. 



(3.4) 



In particular, the imaginary part of the left member vanishes, which means 
that the mentioned parameters are subject to one real constraint. Taking 
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the logarithmic time derivative of (3.4) and using (2.11) give an evolution 
equation for b(t): 



d_ 

di 



/d m d n 

k=l j=l 



We now proceed to evaluate the Poisson integral P((,t) in (3.1). This 
can be done by a simple residue calculus in (2.7), using that \g((,t)\ 2 = 
<7(C>0#*(C>0 1S a rational function in (. The calculation becomes even more 
transparent if everything is done at an algebraic level, by which it essentially 
reduces to an expansion in partial fractions. 

The rational function q{t)/g((,t)g*((,t) has poles at the zeros of g and 
g*, i.e., at u>i, . . . , u m , u\, . . . , io* m . At infinity it has the behavior (by (2.15)) 

lim - V - =A 00 = I l 6 l 2 n j=1 ^ 601 3 6 ) 

^°°9{(,t)g*((,t) [o if m >n. 

Assuming for simplicity that the roots u>i, . . . , u m are distinct it follows that 



m —r m 



9((,t)g*((,t) ^iui k k 



A k A k ( 



_(-U k 1 - WfcC 

where the coefficients A k = Ak(q,b,u\, . . . , u m , Ci, • • • , Cn) are given by 



(3.7) 



A _ g(t) _ g a^-con.K-o) 

for 1 < fc < m. Note that A fc ^ (when / G £>i ocu (!)). 

Now, P(C,t) is by definition (2.7) that holomorphic function in D whose 
real part has boundary values q{t)/g{C,, t)g*(£,t) and whose imaginary part 
vanishes at the origin. The function (3.7) itself certainly has the right bound- 
ary behaviour on 83, but it is not holomorphic in D. On the other hand, 
the two types of polar parts occurring in (3.7) have the same real parts on 
the boundary: 

Re— ^ = Re on 83. 

Q-uJk 1 - w fc C 

Therefore, without changing the real part on the boundary we can make 

the function (3.7) holomorphic in D by a simple exchange of polar parts. 
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In addition, one can add a purely imaginary constant to account for the 
normalization of P at the origin. The result is that 



2Ai 



P (c>i) = A) + g_^ 



(3.9) 



for a suitable constant Aq. Since P(0,t) = J d/j, is real, the imaginary part 
of A Q is given by 



2A- 

lmA = Im 



(3.10) 



For the reflected kernel P*((,t) = P(l/(,t) we get 

2A 3 ( 

1 



c 



(3.11) 



On the other hand, the boundary condition on <9D satisfied by ReP shows 
that 

P((,t) + P*((,t)= ^ 2qit \, r (3.12) 

identically as rational functions. Therefore we find, on comparing (3.9), 
(3.11), (3.7) with (3.12) and using (3.10), that A is given by 



A = A^ + Y,— ■ 



(3.13) 



This also shows that 



JUL A- " pa 



P(0) = / 



(3.14) 



In particular, if m > n, so that A^ = 0, we have 



A, 



Now we evaluate the residues in (3.1): 



d\x. 



(3.15) 



Res^ = 



(P(0 



g'iCY 

9(0. 
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For the poles of 5 the calculation of the residues is simpler because P{() is 
regular there: we simply have 



Res f=c . 



^ s(C) 



= -UjCjPiCj), 



where rij is the order of the pole Q. 
We summarize: 



Theorem 3.2. In the case of simple zeros of g we have the rational dynamics 
d. 4 2A fc ^ 2(A k + Aj) ^ 2A k 

-|, og6 = X + E ^, (3,7, 

fc=l SJ S 

^^ = (m-n+l)A , (3.18) 

where the Aj are given by (3.6), (3.8), (3.13). 

The last equation follows by letting ( — > 00 in (3.3) and using (2.17), 
(2.18) and (3.9). The equation (3.16) was obtained, in the polynomial case 
(n — 0), already in [19]. 



It is useful to observe that the reflected Poisson integral P*(C, t) = P(l/£, t) 
is nothing else than the corresponding Poisson integral for the exterior do- 
main F \ D: 

P*(C,t) = ^-[ 1 f\j + C , - (CeP\D)- (3-19) 



This follows from (2.7) and (3.12) by taking into account that the jump of 
the Poisson integral when ( crosses the unit circle is exactly the right member 
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in (3.12). The orientation of the unit circle in (3.19) is the opposite of that 
in (2.7). Note that P* > in P \ D with P*(oo) = P(0) = / dfi, and also 
that the right member in (3.12) vanishes at the poles Q of g, hence that 

p*(Ci,t) = -p(d,t). 

Next we recast parts of Theorem 3.2 more directly in terms of the exterior 
Poisson integral: 

Theorem 3.3. In the case of simple zeros of g we have 

| iogc fc = p*k, t) - ^(i + £ ~ £ r^r-^ (3 - 20) 

|logG = P*(C„t). (3.21) 



For later use we recall (see (3.8)) that 

2A k 2q 



u k u k g'(u k ,t)g*(u k ,t)' 

Proof. Set 

R(U)- 2q{t) 



(3.22) 



g(C,t)g*((,ty 

so that P + P* = R by (3.12). Returning to (3.1) we then have (suppressing 
t from notation) 

u k = Res[(CP*(C) - CR(C))(log<?(C))'], 

where we have to compute the right member. The residue of the first term 
is immediate: 

Res CP*(C)(log»(C))' = w fc P*K)- 

C=^fe 

For the second term, there is a slight difficulty due to the fact that both 
factors in the product R(()(\ogg(())' have poles at ( = u k . However, we can 
use that 

(\ogRy + (\og g y + (\ogg*y = o 

by the definition of R, hence that 

CP(C)(iog<?(C))' = -CP(C)(iogP(C))' - Cfl(C)0og<7*(C))' 
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= -(R , (o+(- 1 Rmogg(oy* 
= -(Ci?(c))'+i?(c)(i + r 1 (iog^(or) 

Here we used (2.18) in the last step. 

In the last expression, the first term has no residues since it is a pure 
derivative, and in the second term only R(() has a pole (a simple one) at ui^. 
The residue of R(() is 

ResP(C) = f =2A fc , 

hence we get 

Res ci2(c)0ogff(c))' = 2^(1 + E r^77 - E r^)- 

Now (3.20) follows. The other equation, (3.21), is immediate. 

□ 



It is easy to check that (3.21) is valid also in the case of multiple poles (and 
zeros). An immediate consequence is that, as has previously been observed 
in for example [20], [4], [6], poles always move away from the origin. More 
precisely we have 

Corollary 3.1. For each 1 < j < n we have, when q(i) > 0, 

|log|0(t)| = ReP*(0,t)>0 
along with the estimates 

10(01 -i < <* logical \CM+1 ( , ^ 

10(01 + 1 - dlog fll (t) - 10(01 - 1' 1 ' ' 

Proof. It remains only to prove (3.23), and since 

d\og\Q(t)\ _ dhg\CM/dt _ ReP*(Q,t) _ ReP(C*,t) 
dlogai(t) dloga^/dt P(0) P(0) 

this follows easily from the ordinary Harnack inequalities for the positive 
harmonic function ReP(£, t) in D, applied at the point ( = Q. 

□ 
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Because of the middle term in the right member of (3.16), or the second 
term in (3.20), the zeros show up a more complicated behavior than the 
poles, and in particular they do not always move away from the origin (see 
next section). However, when m > n they collectively move out, in fact even 
faster than the poles. This follows by taking the real part of (3.5) and using 
(3.13), (3.14), (3.18) and the definition (3.6) of A^: 

Proposition 3.2. 

d A d A (2(fdfi-A 00 ) if m = n, 

- ^\ogu k {t)-- ^\ogQit) = \ ( .owj - f ^ 

dt ^— ^ dt y{m — n + 2) j dji if m > 

In particular, if m > n, 



71. 



n 



^$>gMOI>^X>g|C;WI, 

k=l j=l 

d m d n 

— ar s u k {t) = j t Y^ ar § o (*) • 

k=l j=l 

The case m = n is indeed a little exceptional in that f((,t) in this 
case has a simple pole at infinity which may be lost at one moment of time 
(see discussion after (2.15)), causing m to temporarily drop below n. What 
happens in (2.15) then is that (at least) one root Ukif) rushes to infinity while 
b{t) goes to zero. Then also goes to infinity. However, after the event 
the root comes back again. 

Example 3.2. The following example of an off-center injection/suction of a 
disk is similar to examples which have been given by Y. Hohlov, S. Howison, 
S. Richardson and others, see [7], [18]. We consider a Hele-Shaw evolution 
with qit) > such that at one particular instant, which we choose to be 
t = 2, the fluid domain is Q{2) = 0(1,2). Then for a general t the domain 
Q(t) will satisfy the quadrature identity 

/ hdm = 2n(Q(t) - Q(2))/i(0) + 4tt/i(1) (3.24) 

JQ{t) 

for functions h harmonic and integrable in the domain. The mapping function 
at t = 2 is 

/(C,2) = 2^C' 
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while it for t ^ 2 can be taken to be of the form 

,C(C -«(*)) 



f(C,t) = b(ty 



for suitable real- valued functions a(t) and b(t). Here we have left the exact 
choice of q(t) open in order to allow the time coordinate simply to be the loca- 
tion of the pole on the positive real axis, which by Corollary 3.1 is consistent 
with the assumption q(t) > 0. The range of t will then be 1 < t < oo. 
The derivative of / is 

e-2tC + ta{t) 
g(C,t) = b(t) ^_ t ) 2 ' 

which has two zeros u>i(t), u> 2 (t) satisfying 

'uw* = ta(t), 

UJ\ + CO2 = It. 

Using the fixed data, contained in the last term in (3.24), one gets the equa- 
tions 

6(1 - at) 
J[ ' j t(l - 1 2 ) ' 

r, ,n b 2 {l-2t 2 + at 3 ){a-t) A 

which can be solved for a, b, giving 



„ m _l , (t 2 -l)(t 2 + l+y/(^- 1)2 + 16) 
a[t) ~ t + 2t(f - 4) 

t 



Kt) = -(t 2 + i-^/(t 2 -iy + i6). 

Now the behavior of the roots can be read off from (3.25): for 1 < t < 2 we 
have U1U2 = at < 0, tending to —00 as t — > 2, while u\ + 002 = It remains 
finite. Also, b(t) — > as t — > 2. Thus the two roots go to infinity along the 
real axis (in different directions). 

For 2 < t < 00, U1U2 = at > 0, which means that the roots come back 
along the imaginary axis as t increases from 2. Finally, as t — > 00, they turn 
back to infinity the same way. 
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4 Global locally univalent solutions are uni- 
valent 



In this section we show that any global locally univalent solution of the 
Polubarinova-Galin equation actually has to be univalent all the time. It is 
intuitively clear that it has to be so. In fact, if univalence breaks down while 
local univalence still holds, then this means that two different parts of /(D, t) 
start to overlap, hence a hole in the part of C covered by /(IB, t) is created. 
With the solution being global this hole eventually has to be filled in, but it 
is easy to realize that when this occurs a branch point for the covering map 
from the multiply sheeted /(IB, t) to C is created. In other words, a zero of 
/' moves into D and local univalence is lost. 

A rigorous proof of the above statement can be based on a result from 
Lowner theory, namely Theorem 6.2 in [15]. To prepare for that, recall that 
in terms of the Taylor expansion (2.10) of f((, t) the zeroth order moment is 
given by 

oo 

M (f)=X>k-(f)| 2 . 

J'=l 

Let 

oo 

*b(0 = Z>;(*)| 2 (4-1) 
i=2 

denote what remains after the first term. Then 

N (t) = M (t) - ai (t) 2 = M o (0) + 2Q(t) - ai (t) 2 

(recall (2.5)). 

In [10] O. Kuznetsova showed that N (t) is a decreasing function of t, 
hence that ai(t) 2 — 2Q(t) is increasing and that 

< N (t) < iVo(O), (4.2) 

aiitf > ai (0) 2 + 2Q(t). (4.3) 

It also follows that 

= = = , ; / 1 as t ->■ oo, 4.4 
y/M {0) + 2Q{t) ^AUt) 
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and that 



|,„ (/) | < for k > 2. (4.5) 



Theorem 4.1. Le£ f(C,t) G 0\ ocu (JD>) be a global (0 < t < oo) locally uni- 
valent solution of the Polubarinova-Galin equation (2.1) driven by injection 
with rate q(t) > such that Q(t) — > oo as t — > oo. Then f((,t) G (9 univ (D) 
for all < t < oo. 

Proof. The choice of time parameter is irrelevant as long as Q(t) — > oo as 
£ — >■ oo. Therefore, by (4.4), we may choose t so that 



ai(t) = e*. 

Then f(C,t) solves the Lowner-Kufarev equation (2.6) with P(C,t) being a 
Caratheodory function, i.e., satisfying 

ReP(C,t)>0 in©, P(0,t) = l. 

This is the setting in [15]. 

Now, for (say) \(\ < | and i> we have, using (4.2), 

oo 

l/(C,0l<e' + £M*)|2-'" 
i=2 

(oo \ V 2 / oo \ V 2 

X>'M*)I 2 J (Ej 2 ^' 

< e* + v/iVoW < e* + ^(0). 

Thus |/(C, *)l < C e * (0 < t < oo) in a neighborhood of the origin, which 
means that the assumptions in Theorem 6.2 of [15] are satisfied. The con- 
clusion is that /(C) is univalent for all t > 0. In our case it follows that 
/(C) is actually univalent in the closed unit disk, because by assumption it 
is locally univalent there, and if for two different points (1,(2 £ <9D we had 
/(Ci, = /(C2, £) for some t then univalence in the open disk D would be lost 
in the next instance (the boundary is really propagating with positive speed 
under our assumptions). □ 

Remark 4.1. Under the slightly weaker assumption that /(C, t) G (9 norm (D)n 
Ciocu(D) (i-e., local univalence only in the open disk) the theorem holds with 
the conclusion correspondingly changed to f(C,t) G O un 
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5 Asymptotic questions in the polynomial case 



In this section we study asymptotic behaviour of the roots of g — f for 
locally univalent polynomial solutions of (2.1) with q{t) > and Q(t) — > oo 
as t — ?■ oo. Thus n = and we write 



We assume a m+ i(0) ^ 0, by which a m+1 {t) ^ for all t, since M m = 
a m +i{t)ai{t) m+l (see (5.3) below) is preserved. 

5.1 Collision of roots 

In the rational case we have seen (Example 3.2) that the zeros of g may 
occasionally (when m = n) move to infinity and come back again, and we 
have seen that the poles always move out. In the polynomial case (n = 
0) the zeros cannot move to infinity in finite time, as can be seen from 
Proposition 3.2 for example. Below we give some examples showing that 
zeros can still move towards the origin and that they can collide, even though 
they in principle repel each other (as can be seen from (3.16) for example). 

Example 5.1. When m = 2, n = and the mapping function / has real 
coefficients the solution of the Polubarinova-Galin equation can be made fully 
explicit (see Example 5.2). The roots uj\ and U2 will either be a complex 
conjugate pair, or both will be real, with occasional collisions allowed. We 
shall first make some immediate conclusions from the dynamical equations 
(3.16) in the case that the roots oui, 0J2 are real. 

Assume for example that both roots are positive, say 



m+l 




(5.1) 



1 < U>i < U>2- 



By (3.8) the coefficients A k are then given by 




(uj 2 - wi)(l - w|)(l - wxwa) ' 



(Ui - W 2 )(l - wf)(l - U-1U2) ' 
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Now (3.16) together with (3.13) gives 



\b\ 2 d . . 0J1U2 

logwi = a{ui,u} 2 ) + P{uji,uj2) 



where 



q dt U2 — u>i 

\b\ 2 d UJ1U2 - 3 

— \ogu 2 = a{ui,u 2 ) - p{u u u 2 ) , 

q at co 2 — oji 

2(1 + U1U2) 



K-i)K 2 _i)(^ 2 _i)' 

(wi + w 2 ) 



(a;? - l)(a;| - l)(a; 1 a; 2 - 1) ' 

Note that a(wi, w 2 ) and f3(u\, u> 2 ) are positive and symmetric. The remaining 
factor ^^J 3 has a singularity when ui\ = 0J2 and changes sign when U1U2 = 3. 
Hence we can draw the following conclusions. 

• The product U1U2 always increases in time (also clear from Proposi- 
tion 3.2). 

• The ratio u> 2 /u>i decreases if Uiuj 2 < 3, increases if Wiw 2 > 3 

• If < 3 and U2 — oo± is sufficiently small, then 6j\ < and 002 > 0. 
Hence no collisions occur for Wiw 2 < 3. 

• If = 3 then u>i = u 2 > 0. 

• If U1U2 > 3 and 0J2 — u>i is sufficiently small, then ui > and C02 < 0. 
The condition 0J1U2 > 3 is preserved in time and eventually leads to a 
collision (even if 0J2 — oo± is not small to start with). After the collision 
the roots cease to be real and continue as a complex conjugate pair. 

In the case of complex roots (but still with / having real coefficients) the 
scenario may be quite interesting, as the following example due to C. Hunt- 
ingford [9] shows. 

Example 5.2. When expressed in terms of the moments Mi, M 2 and with 
the time parameter chosen so that ai(t) = e* (q(t) > defined accordingly), 



22 



the solution of the Polubarinova-Galin equation (2.1) in the case m — 2, 
n = with / having real coefficients is explicitly 




By taking 



one gets an interesting example which has been investigated by C. Hunting- 
ford [9]. We let the solution start at the first possible instant of time, namely 
t — to — | log | < (for smaller values of t, f is not locally univalent). At 
that moment the image domain Q(t) = f(3,t) has two cusps on the bound- 
ary (the two zeros of g form a complex conjugate pair on the unit circle). As 
t increases these cusps resolve, the zeros of g move out from the unit circle 
and collide on the real axis at some moment t — t± < 0. After that one of the 
zeros, say u>\(t), moves back to the unit circle and reaches it again at time 
t — 0. Thus a new cusp (which will be a 5/2-power cusp) develops on the 
boundary. However also this cusp resolves, the root oj\(t) moves away from 
the unit circle along the real axis, captures and collides with the other root 
0J2{t), after which they leave the real axis and finally move towards infinity 
in the asymptotic directions of the positive and negative imaginary axes. All 
of this can be seen by examining [9] carefully. The final asymptotics follows 
from Theorem 5.1 below. 

The solution is global (t < t < oo) and is all the time univalent (oth- 
erwise it would have broken down in finite time, by Theorem 4.1), or the 
remark following it. 

The above example shows that it is possible for a single zero to reach the 
unit circle in the injection case. The next example will show that this can 
not occur if the other zeros are sufficiently separated and far away. 

Example 5.3. Let n = 0, m > 2 with / allowed to have complex coefficients. 
We assume that one of the roots is close to the unit circle and the others far 
away and well separated from each other, say 

1 < ui < 1 + e, e > small, 

|w fc - u)j\ > M (1 < k, j < m, k ^ j, M > 1 large). 
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We shall use equation (3.16) to investigate short time root behavior. 
First observe by (3.8) that 



A 



k 



u k |6| 2 (i - M 2 ) n^ fe K-^)K-^)' 

hence that, among these quantities, ^ is the dominating one under the 
present assumptions. Therefore, by (3.16) and (3.13) (with = 0) we have 

-— log cui (t)\ = ReA + Re — - + Re > ^ 

at uji u>i — ujj 

3Ai 2A 1 ^ ui , 
« Re h Re > J 

3Ai 3g „ 1 

« Re — - = ^tt — ^ — • Re 



ito-m 2 ) nr =2 (^-^-)K-^) 

3g 



i^i 2 (! — i^ii 2 ) n7^2 i^i — 



<0, 



L 3= 

which shows that the minimum root u>i(t) moves away from the origin to 
prevent the strong solution from blowing up in short time. 

Similarly, in the case of suction, under the same assumptions the mini- 
mum root u>! moves towards the origin initially and will make the solution 
blow up. 



5.2 Long time behaviour of roots and coefficients 

We now assume that the locally univalent polynomial solution (5.1) is global 
in time. By Theorem 4.1, f((,t) is then actually univalent all the time. We 
shall describe the asymptotic behavior of roots of g((,t). It will be shown 
(Theorem 5.1) that these roots all move away from the origin as time t 
is large, even though some of them may move towards the origin initially. 
One naturally expects that what makes the initial root behavior of g(C,t) 
unpredictable is that the distribution of zeros of g((, 0) is not always uniform. 
Therefore, it seems reasonable that, by controlling the distribution of these 
initial roots, we can guarantee that all roots always move out. In fact, this 
will be demonstrated in Theorem 5.2. 
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Recall first Richardson's formula [17] for the harmonic moments M k (see 
(2-4)): 

M k = ■ ■ ■ a i k+ i a h+-+i k+ i- ( 5 - 2 ) 

(u,— ,ik+i) 

Here k — 0, 1, 2, . . . and the summation runs over all (/c+l)-tuples (ii, • • • , ik+i) 
of integers with ij > 1, with the convention that = for j > m + 1 in the 
present polynomial case (5.1). See [17] and [10] for more details. Note that 
the final nonzero moment has a quite simple expression: 

M m = < +1 <wI. (5.3) 

Lemma 5.1 below will show that asymptotically, as t — > oo, the formula (5.3) 
will almost be true for all moments. 

Besides the zeros Uk(t) of g(C,t), which are what we want to study, we 
introduce two other sets of zeros, for comparison: 

• Let 

tik(t) = uJkifya-tit)-^ (5.4) 

m+2 

be corresponding rescaled zeros, namely the zeros of g((, i) = g(a 1 m (, t). 

• Let 

1 

Wfc = 



y/-(m + l)M m 

be the zeros of the polynomial a± + (m + l)a™ +2 a m+ iC m consisting of 
the first and last term in g((,t). 

Note that the Uk do not depend on time. On the other hand, the 
rescaled zeros Uk(t) are normalized in such a way that they asymptotically 
stabilize, and approach the Co k (Theorem 5.1). Note also that the product 
cDi(t) • • -Cjmit) = — , TTvff is a conserved quantity. 

The following lemma is a slightly sharpened version of estimates first 
obtained in [11]. 

Lemma 5.1. Given m > 1 there exist universal polynomials Fj(x±, . . . , x m+2 ), 
1 < j < m + 1, in m + 2 real variables such that 

\ aj (t) ai (ty -M~i\< ^-M^y V^o(bT> |M|, • • • , \M m \) (5.5) 
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(t > 0, 1 < j < m + 1 ) whenever f((,t) given by (5.1) is a locally univalent 
solution of (2.1). The polynomials Fj are increasing functions in all their 
variables. 

In particular we have estimates of the kind 

IMWI < (5-6) 

Mt)\<-^g. 2<j<m+l, (57) 



for suitable constants Cj = Cj(^oj, v^o(O), \Mi\, |M 2 |, • • • , |M m |). 

Proof. Using (5.2) we first prove by induction on decreasing j that there exist 
polynomials Gj(x 1 , . . . , x m+2 ), l<j<m+l, inm + 2 variables, increasing 
in all variables, such that 

\a 3 {t) ai {ty -M~i\< -L^-^, |Afi|, • • • , |M m |). (5.8) 

We start by setting G m+ \ = 0, which makes (5.8) hold for j — m + 1. Now 
we show (5.8) for j = k assuming it holds for all j > k + 1. The arguments 
will work for any 1 < k < m. Using (5.2) and suppressing dependence on t 
from notation, we have 

M k _i — a k a\ = ^ iia ix ■ ■ ■ a ik ■ a il+ ... +ik 

(h,- ,ifc)^(l,...,l) 



E 



l t a n ' " a ik n h+-+ik 

ii+-+ifc ' a n+-+ik a l 
(h,- ,i k )^(l,...,l) "1 



1 v-^ hah ■ ■ ■ cti k il+ ... +ik 

2 _ii+-+ifc-2 ' "w+-+*fc U l 
1 (n,...,i fc )^(l,...,l) fl i 

We shall estimate the terms in the above sum. Since (ii, ■ ■ ■ ,i k ) ^ (1, . . . , 1), 
at least one of the ij is > 2, hence it follows from (4.5), (4.3) that the first 
factors can be estimated as 



'h'^Qj^-^ • • • 



iiH hjfe-2 



0<o,/3<m+l ±v ' 
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Moreover ii + ■ — h i k > k + 1, so by the induction hypothesis we have 



I h-\ hifc I 

\a il+ ... +ik a l — ivi il+ ... +ik _i\ 

1 



< G il+ ... +ik (— -, v/iVoM, \M,\, • • • , \M m \)- 



ai(0) ai(t) 2 
in particular 

In- - n il+ '" +ik \ 

|"iiH h« fc u l I 

1 1 m 

< -, v/iVO), IMxI, • • • , \M m \)— — + Y, 

Therefore we can estimate also the second factors in the above expression for 
a k a^ — Mfc_x. From this we easily deduce (5.8) for j = k knowing that it is 
true for j > k + 1. This completes the induction step and hence proves (5.8) 
for 1 < j < m. 

Now, already (5.8) shows that the estimate (5.7) holds. Since this esti- 
mate improves (4.5) by a factor at least ai(t) 2 in the denominator, we can 
'bootstrap' the previous argument: using the new estimate in (5.9) makes the 
induction process work with the factor l/ai(t) 2 in (5.8) replaced by l/ai(t) 4 . 
Thus (5.5) follows, and since (5.6) is just the special case j = 1 the lemma 
is proved. 

□ 

Theorem 5.1. With notations as above, assume that /(C,t) is a global poly- 
nomial solution of (2.1). Then 

\uJk{t) — 0Jk\ — > as t — )• oo 

if the roots Uk(t) and Cj k are ordered appropriately. Furthermore, all roots 
eventually move away from the origin as time is large enough. 

Proof. The monic polynomial vanishing at the rescaled roots (5.4) is 

n«-%)=E(^) jd^-x 1 

j=i 3=1 



(m + l)a m+ i 



y, 1 j ajaj -±(m+l-j) .j.! 

f^m + l a m+1 a? +1 1 



\ , v 3 ajai n -™ (m+1 - j) n- 1 4. r 



(m+l)M m ^m + lM, 



27 



Due to Lemma 5.1 and because the exponents —-^{m + 1 — j) are strictly 
negative the coefficients of the middle terms tend to zero: 



J 



QjjQj-y 



m + 1 M„ 



< 



J 



-Am+l-j) 



m + 1 a\M n 



+ 



3 



M. 



J'- 1 . 



-(m+l-j) 



m + 1 M„ 



->■ as t ->• oo (2 < j < m). (5.10) 



From this it follows that | — ci}fc| as t — >■ oo. 

It also follows that the roots move away from the origin as time is large, 
because the speed of the roots only depends on the position of the roots, 
and it is clear that for a symmetric configuration of roots, like Uk, the speed 
points radially away from the origin. A slightly more precise argument can 
be based on (3.20), (3.22), by which 



^log|w fc | = ReP*(w fe ) 



Re 



2q 



:i + 



m 



1 



^1- 



Here the first term is always positive, while the subtracted term is nega- 
tive whenever the roots are sufficiently far away and close to symmetrical 
configuration. In fact, in this case Ukg'{o0k) ~ —mai, g*{oJk) ~ a \ an d 



i + E 



i 



j=i i- 



1. (This argument will be made more precise in the proof 



of Theorem 5.2 below.) Hence 4 log \ujk\ > for large t, as claimed. 



□ 



5.3 Assumptions only on initial data 

In this subsection we do not assume a priori that f(C,t) is global, we only 
make assumption on the initial data, and the solution being global will be 
part of the conclusion. 

Theorem 5.2. In terms of any initial function f((,0) = Y^f=i a j(0)C' 7 ; ^ 
M > be a common upper bound for |Mi|, . . . , \M m \ and ^/JVo(O). Then, for 
any e > 0, there exists a number B = B(e,M,m) > such that whenever 
ai(0) > B the polynomial solution f(C,t) starting with f((,0) is global in 
time and the assertions (i)-(iii) below hold for 1 < k < m. 
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(i) For a suitable ordering of the ui k {t) and Cj k we have 

\CJk{t) — &k\ < e i < t < oo. 
In particular the roots u k never collide (with e > sufficiently small). 

(n) 

\oJk(t) — cu k \ — > as t — > oo. 
(Hi) The roots uj k {t) move away from the origin (for all t > 0). 

Proof. First, choose p(M m ,m) > so that the disks 3(6jj,p) (1 < j < m) 
are disjoint and such that, for any point z G 3((jj,p), 



z m + 



[m + 1)M„ 



< 



A(m + 1)\M„ 



■ tan • 



7T 

20' 



(5.11) 



tan t|j will do. 



For example, any < p < \\M m 

Next, as in the proof of Theorem 5.1 we have that for any < e < p 
there exists 5 > such that \u k — Cj k \ < e whenever 



-(m+l-j) 



m + 1 Mr, 



< 8, 2 < j < m. 



(5.12) 



Finally, it follows from Lemma 5.1 that (5.12) indeed holds if just cti(0) is 
large enough. We conclude that for any < e < p, there exists Bx(e, M, m) > 
such that if ai(0) > B 1: then \u k (t) — Co k \ < e and (5.11) holds. Therefore 
(i) holds. By Theorem 5.1, also (ii) holds. 

From (3.20), (3.22) we see that in order to prove that all roots move away 
from the origin, it is sufficient to prove that 



Re 



u k g'(u k ,t)g*(u k ,t) 



1 



j'=i 



UjtU k 



< 0. 



(5.13) 



We first estimate the terms of u k g'(u k ,t) = ^7=1 JO + ^) a j+i u k- F° r 1 ^ 
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j < m — 1 we have 

m / 

1 / 

_ j+i 2 (4- 1 ) 

Hence, for 1 < j < m — 1, 

KW| < [1 + (e + l^l)™" 1 ] (m + ^ a" 2/m 0l (5.14) 

with Fj = ^j(^To), M, • • • , M) as in Lemma 5.1. From (5.14) and (4.3) we 
conclude that there exists B 2 > B\ such that if a±(0) > B 2 , then + 
l)a j+1 ul\ < |aitan^. For j = m, a j+ iu J k = a m+ iuj™ = u^Mjai by (5.4). 
When ai(0) > -B 2 we have 

uj k g'(uj k , t) + mai = ( u k + — — J M m m(m + l)on 

V M m (m + 1)/ 

m— 1 

due to (5.11), and since + l)a,j + iu 3 k \ < |aitan^ we conclude from this 
that | &rgu k (t)g'(uj k , t) — n\ < ^. Finally, we can find B = B(5,M,m) > 
B 2 such that if ai(0) > B, then | aigu k (t)g'(u k , i) — n\, \ &rgg*(u k , t)\ and 
| arg(l + J2J=i ^r)l are an < is- Therefore (5.13), and hence (Hi) in the 
theorem, holds. 

□ 





6 Asymptotics for rational solutions 

In this section we study the asymptotics of poles and Taylor coefficients in 
the rational case. We assume that / G 0i ocu (O) and g — f are of the forms 
(2.20) and (2.15) and that q(t) > with Q(t) -> oo as t ->■ oo. Recall (3.17), 
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(3.21) and Corollary 3.1, which in particular show that the poles always move 
away from the origin. The following theorem gives more precise estimates of 
their locations. 

Theorem 6.1. Assume f((,t) is a global in time locally univalent solution 
of the Polubarinova-Galin equation, with g((,t) rational of the form (2.15). 
Then \Q{t)\ ~ a\(t) for each 1 < j < n. More precisely, with ai(t) as in 
(2.10) 

1 (IO(0)l + ^-2)<^#<-4v(|C,(0)| + -i^ + 2) 



ai (o)^ JK 71 IG(o)| " - ai (t) - Ol (0) vwv 21 10(0)1 

for allt > 0. 

Proof. Setting r = logai(t), £ = £(r) = \(j(t)\ > 1, the Harnack estimates 
(3.23) say that 

m^l < g < ttL±H 

f + 1 - dr ~ £ - 1 

On integrating the differential equations corresponding to the equality cases 
one easily obtains the inequalities 

^ (ICj(0) | + ^- 2) ^ ( i at) | + ^- 2) 

and 

^^' + W + 2)s 4 (l&(0),+ KM + 2) ' 

Since 



^■ (lc ' (t)l+ w- 2)£ ^-^ (l0(t)l+ w +2) 

the desired estimates follow. 

Next, turning to coefficients we shall write / and g as follows: 



□ 



/(c,o = $>>g(i - 7 )-) + ^ ;;;; = £%wc j , (e.i) 
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g(c,t) = ±ry . ' = Yti + i)oi+i(*)c j - (6.2) 

Introduce also the Taylor coefficients 5j of the pure rational part of f(C,t) 
by 

Tn 1 -. -E^oc (6-3) 

The above expansions are to be compared with (2.15) and (2.20). In 
particular, £ denotes the number of different (finite) poles of g, rij denotes 
the order of the pole at Q (as in (2.15)) and 



3=1 3=1 



is the total order of the finite poles of g. The numbers Cj are the residues of 
g(C)dC at the points Q, and some or all of them may be zero. 

At infinity, g(()d( has a pole of order no = m — n + 2 (as a differential), 
hence / has a pole of order m — n + 1 there. This is also what (6.1) gives, 
hence the notations in (6.1) and (6.2) are consistent with those in (2.20) 
and (2.15). Recall also that the coefficients ej do not depend on t. The 
coefficients bj, Cj, dj, bj, Cj satisfy 6 m -^+i 7^ 0, c n _^ 7^ 0, b\ 7^ 0, c 7^ 0, 
b m 7^ 0, c„ 7^ 0, di 7^ 0, &i 7^ 0, c 7^ 0, and we shall normalize them so that 
c = 1, c = 1. 

Lemma 6.1. Assume that f((,t), given by (6.1), is a global solution. Then, 
as t — > oo ; 

\c n -i\ ~ a 1 ^ ^, \cj\ = 0(a^) for 1 < j < n — t — 1, 

|&i| ~ fii ~ ai, |6j|=0(l) for 2 < j < m -£+ 1, 
|c n ] ~ a]~ n , = 0(a7 J ) for 1 < j < n — 1, 

|&o| ~ ai, = 0(1) for 1 < j <m. 

See Subsection 2.1 for the meaning of ~. 

Proof. Since Ej=o c jWC' ? — c n-e(t) YYjZiiC ~ Cj(t)) an d c o = 1 the estimates 
for Cj(t) follow immediately from Theorem 6.1. Similarly for Cj(t). 
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Since by (6.1), (6.3), 

t 

we have 

|%(t)-a,(t)|<^-->0, 

as t — > oo. The coefficients bk are given by 

fc-i 

bk = SiCfe-i + ^2 ®3 c k-j + 5 fc, (6.4) 
i=2 

and since the cij, and hence the aj, j > 2, are bounded (see (4.5) for example) 
the assertions about the bj follow easily. For the bk we have similarly 

fe-i 

h = axCk + ^(j + l)a j+1 c k -j + (k + l)a k+u (6.5) 

3=1 

hence the estimates for these follow in the same way. 

□ 

Theorem 6.2. Assume that f((,t) in (6.1) is a global solution and introduce 
the truncations 

N 

f N ((,t) = J2 a At)( J - 

3=1 

Then the following assertions hold. 

(i) There exist numbers s N with s N — > oo as N — > oo such that, for 
J = 0,1, 

fjP(C,t)-f®(C,t) =0( ai (t)- s »), (6.6) 



sup 



as t — > oo. The same is true for any j > 0, with sn then depending on 
j- 



(ii) For each k > 2, 



lim ak(t) ai (t) k = M k -i. (6.7) 

t— >OD 
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(Hi) Assume that, for some number r > 2, Mi = ■ ■ ■ = M r _i = 0, M r ^ 0. 
Then 

lim a s {t)a 1 {t) r+l = 0, 2 < s < r, (6.8) 

and r < m. In case there are no logarithmic singularities (i.e., ej = 0, 
1 < j < £) then we even have r < m — t. 

Proof. Write / on the form 

/K ,o=|>>g(i-|) +6 4±^§ 

oo / £ \ oo 

= - E E tk K + ^ E(^(o fc + piemen 

fc=l \j=l / fc=0 



namely with 

m 

6i 



Then all coefficients in P and Q are (^(a] -1 ) by Lemma 6.1, hence 

sup |P| < — , sup \Q\ < 

O(0,R) a l D(0,R) a l 

for any fixed R > 1 and suitable constants C R . It follows from the above 
that for large N the remainder 

oo 

/(C) - MO = E "A ' 

j=N+l 

is built up by terms of the kind (Q(() k and (P(()Q(() k with also k large. 
Considering only t so large that < \ (for example) it follows that 

sup |/ - f N \ < % (6.9) 
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where s can be made arbitrarily large by choosing iV sufficiently large. Now 
(6.6) follows. 

(ii) Recall the expression (2.4) for the moments M k = Mk(f) in terms of 
/. The truncations f'jy of / similarly define moments: 

M k (f N ) = ^-[ ti(C,t)f N (c,t)fM,t)dC 

Z7n Jan 

— ^ ] HQ-ii ' ■ ■ a i k+1 a h-\ Nfc+i- 

It is clear from (2.4) that each moment M k (f) is a Lipschitz continuous 
function of / if / is measured by the norm supp(|/| + |/'|), or any norm 
su PD(o,i?) I /I) R > 1- Therefore we have, by choosing N = N k large enough in 
(6.6), with j — 0,1 and on using (6.9), that for any exponent p > we have 

|M fe (/)-M fe (/ JV )|<C"|/-/iv|<^ 

for some constants C" and C" (which depend on / and p). From this we get 
(with M k = M k (f)) 

at +1 a^I - M k 
= M k (f N ) -M k - hoi! ■ ■ ■ a ik+1 a h+ ... +ik+1 

k+2<Y / ij<N 

= (~p)~ Yl ii a h ■ ■ ■ a i k +i a h+-+ik+i ( 6 - 10 ) 

1 k+2<J2ij<N 

= °{-p) ~ ( k + 2)a^a 2 a^- ^ ha h ■ ■ ■ a ik+1 a il+ ... +ik+1 . 

1 fc+3<£ij<JV 

Now, to prove (6.7) we shall prove by induction that for every s > 2, the 
two assertions 

lim a^tyajt) = M s _ 1: (6.11) 

t— >oo 

\ a M = 0(^-y g ), J>S (6.12) 

hold. First, on using the fact that the |%(i)|, j > 2, are uniformly bounded 
we deduce from (6.10) that 

\a k , + W| < \M k \ + 0(\) + 0(4), k > 1, 
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and hence that |afc + i| = O(^j-), k > 1. Now repeating (6.10) with the new 
estimate \a,j\ = O(j^), j > 2, we obtain 

K + W| < |M fc | + O(^) + O^ 2 ), fc > 1, 

and hence |afc+i| = O(t-), > 1. 

To start the induction process, take k = 1 and use the last estimate 
\ aj \ = O(^), j > 2, in (6.10). This gives 

a^_ Ml = (l) + 0(4) 

and hence lim^oo afa2 = M 1 . Therefore (6.11) and (6.12) hold for s = 2. 

Now take Sq > 2 and assume that (6.11), (6.12) hold for all s < s . Then 
we shall prove (6.11), (6.12) for s — s + 1. Thus we may in (6.10) use the 
fact that for any s < s we have \aj\ = O(^), j > s. This gives 

\a1 +l a k+1 \ < \M k \ + 0(\) + Oiat 2 - 30 ), k > s , 

a i 

and hence |a fc+ i| = 0(-^ TT ) for k + 1 > s + 1. It follows that (6.12) holds 
for s — so + 1. 

Using, in (6.10) with k = Sq, that \aj\ = 0( S p+i ) for j > s + 1 (just 

proved) and \a,j\ — 0{\) for 2 < j < s (induction hypothesis), we obtain 

a i 

and hence lim^oo a^ 0+1 a so+ i = M so . Therefore also (6.11) holds for s = s +l. 
Thus (u) in the theorem is proved. 

(iii) Assuming now M x = ■ ■ ■ = M r _ x = 0, (6.10) with 1 < k < r - 1 
gives 

^ = °(^Tk+r)- Yl ll ~ k+i k+1 a h+-+ik+i ( 6 - 13 ) 
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with p arbitrarily large. Using (6.7) and choosing p > 4, the above gives 




With k = r — 1 this gives (6.8) for s = r. It also follows for s = r — 2. 
Assume now, as an induction hypothesis, that (6.8) holds for s < s < r. 
If s > 3 we then let k = s — 2 in (6.13), and use (6.7) and the induction 
hypothesis. This gives (6.8) for s — so — 1. Hence (6.8) is proved. 

The identity (6.5) holds for 1 < k < m, but also for k > m with the 
convention that b k = when k > m, and Cj = when j > n. Choosing 
k — m + 1 then gives 



If r > m + 1 then (6.8) and Lemma 6.1 show that the right member is 
o(a^ r+1 ^ ■ 0(a^ 1 )) = o(a^ r+2 ' ) ), while (6.7) shows that the left member is 
exactly 0(a^ r+2 ^), if M r ^ 0. This contradiction shows that r < m. 

Similarly, using (6.4) one obtains r < m— t in case there are no logarithmic 
terms in /. An alternative way of proving these upper bounds for r is given 
in the remark below. 



Remark 6.1. The fact that for rational g, the vanishing of a sufficiently long 
sequence Mi, M 2 , . . . , M^ of moments implies the vanishing of all Mj, j > 1, 
can also be deduced from the relevant quadrature identity, like (2.14), holding 
for the image domain. Assume for example that g(()d( has no residues, so 
that there are no line integrals in (2.14). Then choosing h(() = -^7 with 

z e C \ Q in (2.14) gives the identity 



rn 



(m + 2)a m+2 



+ l)a j+ i5 m+w . 



□ 




zk+l UtA z ~^ k ET=o +1 c^ 



M k _* n ^ ^ - T^Bjzi 



for suitable Bj, Cj. Recall that J2j=o( n j — I) = m — £ + 1, see at (2.14) and 
(2.20). If now Mi = . . . M r ^ = 0, M r ^ 0, then we get 



(M + — + — — + ... ){C m - e+ i + —— + ■■■ + —zjtj) 
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Br, 



Bm-e-i 



+ 



Bn 



vvn—l ' 



Z Z'' 

Here it is easy to see from the general structure (2.15), (2.20) that Mq ^ 0, 
Cm-e+i 7^ 0, B m -t ^ 0. On the other hand, Cq — C± — • • • — C no _i = 
because the rational function above has a pole of order (exactly) no — 1 at 
z = 0. Thus, in view of the fact that m — £ + 1 — (n — 1) = n — £ the left 
hand side of the above equation actually is 



(m + Mr + Mr+1 + vr + a 

\ M + — + y r+l + • • • ){^m-l+l + 



m—£ 



c 



M U m -i + i H h 



+ 



M C, 



C m -f 



y n-e 



%m— £+1— (no— 1) 

+ . . . 



r+1 • 



2 Z" -1 Z' 2' 

Comparing with the right hand and using that n — £ < m — £ it follows that 
r cannot be larger than m — £, which is also what Theorem 6.2 tells. 

In a similar way one proves that r < m in the presence of logarithmic 
terms in /. 

From Theorem 6.2, we obtain the following result: 

Corollary 6.2. Let f((,t) be a global solution of the form (6.1) and assume 
Mi = • • • = M r _! = 0, M r ^ 0. T/ien, 



lim sup 



/(CO - >/2Q(t) + Mo(o)c] (v^gft))*" 1 - M r C r+1 



0. 



Proof. By (4.1), (4.2), (4.3) and Theorem 6.2, 

-N (t) 



(t) - y/2Q(t) + M (0) 



Ol 



1 



ai(t) + V2Q(t) + M (0) a 



2r+3- 
1 



Applying once more (4.2), (4.3) and Theorem 6.2 the corollary follows. □ 
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